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Spectra of massive QCD Dirac Operators from Random Matrix Theory: 
all three chiral symmetry breaking patterns* 
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The microscopic spectral eigenvalue correlations of QCD Dirac operators in the presence of dynamical fermions 
are calculated within the framework of Random Matrix Theory (RMT). Our approach treats the low-energy 
correlation functions of all three chiral symmetry breaking patterns (labeled by the Dyson index f3 = 1, 2 and 4) 
on the same footing, offering a unifying description of massive QCD Dirac spectra. RMT universality is explicitly 
proven for all three symmetry classes and the results are compared to the available lattice data for f3 = 4. 



1. Introduction 

Random Matrix Theory (RMT) has turned out 
to be a very fruitful tool in studying the phe- 
nomenon of chiral symmetry breaking in low- 
energy QCD 0. First proposed as a purely phe- 
nomenological approach, it has recently been put 
onto firm field theoretic grounds after the analytic 
RMT predictions have been reproduced within 
the framework of finite-volume partition func- 
tions and partially quenched chiral perturbation 
theory using supersymmetry [0 and replica || 
techniques. 

Similarly to previously studied sum rules @, 
RMT solutions for spectral statistics serve as a 
more detailed test of QCD lattice data, in par- 
ticular for a given sector of topological charge v 
(which counts the number of zero modes of the 
QCD Dirac operator). This has become possible 
due to recent developments in lattice gauge the- 
ory. Namely, the Ginsparg- Wilson relation was 
shown to provide an exact chiral symmetry on 
the lattice together with a well defined topologi- 
cal charge v || (and these proceedings). The clas- 
sification of the three different xSB patterns ac- 
cording to gauge group and representation IJ has 
been confronted to lattice data and good agree- 
ment has been found for the spectral density and 
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distribution of the smallest Dirac eigenvalue for 
massless lattice data of all xSB patterns and for 
different values of v (see e.g. Q, Fig. 2). 

It has to be mentioned that the direct field the- 
oretic calculation of spectral correlators is much 
more cumbersome as compared to RMT. The for- 
mer approach has not led to any explicit analytic 
results beyond the spectral density due to enor- 
mous increase of dimensionality of auxiliary (su- 
permatrix) fields involved. So far, only spectral 
densities for massless flavors (in all three xSB 
patterns) and for one single massive flavor (in 
SU(N C > 3) in the fundamental representation) 
have been derived In contrast, the classic 
RMT technique is free of the above technical com- 
plications and allows computing the higher order 
correlation functions with the same ease (see Ref. 
|| for comparison with lattice data). This ad- 
vantage of a RMT description becomes even more 
significant for dynamical fermions. The correla- 
tion functions with an arbitrary number of mas- 
sive flavors Nf have been calculated for funda- 
mental SU{N C >3) §,0 09 = 2) and recently for 
the two remaining xSB patterns JTT| , |T2[ ] {(3— 1, 4). 

In the present communication, we report on our 
results |llj relevant for gauge groups SU (2) in the 
fundamental ((3 = 1) and SU(N C ) in the adjoint 
representation (/3 = 4). Our predictions are com- 
pared to SU(2) lattice data |r|] for dynamical 
staggered fermions with 4 degenerate flavors. 
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2. RMT results for massive flavors 

Let us briefly recall the connection between 
RMT and low-energy QCD. The Dirac oper- 
ator spectrum at the origin is related to the 
chiral condensate S, the order parameter of 
%SB, through the Banks-Casher relation |l4|] 
£ = limy^oo Trp(Q)/V , where V is the Euclidean 
space-time volume. Here, the spectral density 
of the Dirac operator is given by the average 
p(X) = (^2 k S(X — Afc)} over all gauge field config- 
urations and the A& are the Dirac operator eigen- 
values. In the limit A -1 <C V <C m" 1 H , where m„ 
is the pion mass and A is the scale of the lightest 
non-Goldstone particle, the QCD partition func- 
tion is dominated by zero momentum modes of 
the Goldstone fields and hence collapses into a 
simple group integral As a result, the parti- 
tion function only depends on the global symme- 
tries of the QCD Dirac operator and contains just 
the rescaled quark masses pf — mfVYi as parame- 
ters. In a sector with fixed topological charge v, it 
coincides with the corresponding RMT partition 
function once the space-time volume V is iden- 
tified with the size n of the corresponding ran- 
dom matrix. Analogously, all correlation func- 
tions can also be computed from RMT provided 
the matrix eigenvalues are appropriately rescaled, 

= AfeFS = \knirp(0), where V — > oo or n — > oo 
is taken; p(0) denotes the RMT spectral density. 
This provides us with a parameter free prediction. 

The joint probability density function of chiral 
RMT associated with Nf massive quarks in the 
sector of topological charge v is defined as 
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Here, (3 = 1, 2 and 4 labels the symmetry of the 
matrix ensemble to be orthogonal (/3 = 1), uni- 
tary (/3 = 2) or symplectic {(3 = 4) in correspon- 
dence with the three xSB patterns [||. The par- 
tition function appearing in the normalization is 
obtained by integrating over all eigenvalues A&. 
The fc-point correlation function is determined by 



integrating over n — k eigenvalues only: 
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Here, A„ ({A}) = n™<j(-^ — -N?') * s tne Vander- 
monde determinant and the weight function is 
given by 



, y (A)=A^ + f- 1 e- /3V W, 



(3) 



where V(X) is a finite-polynomial confinement 
potential whose form is not fixed apriori. Al- 
though the simplest choice V(X) = S 2 A defining 
the Gaussian ensemble leads to significant mathe- 
matical simplifications it cannot be derived from 
the QCD Lagrangian. It is therefore crucial to 
show that the RMT results for the rescaled fc- 
point correlation functions, Eq. (||), are universal 
and do not depend 3 on this choice for V(X). 

In the following we present a unified way to 
explicitly calculate, and prove, the RMT univer- 
sality of massive spectral correlators for all three 
xSB patterns, (3=1,2 and 4. Our strategy is to 
express the massive spectral correlators in terms 
of the known massless ones; the latter have al- 
ready been shown to be universal Jig] , To pro- 
ceed, we assume that the massive fermions are 
/3-fold degenerate. With the help of the identity 

A„ +JV/ ({A},{-m 2 }) 



A Nf ({-m*}) 



N f 



= A n ({X})lU[(X l+m }) 

i=i/=i 
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associated 



the joint probability density P, 
with the /3-fold degenerate massive fermions of 
total amount (3Nf can be rewritten through the 

massless joint probability density f„+jy with n 
positive {Xi} and Nf negative {— mj} eigenval- 
ues. This leads us to the remarkable identity 
which holds in all generality for finite n: 
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(5) 



3 The only condition is that the macroscopic RMT spectral 
density has to obey p(0) ^ 0. 
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In order to compare with QCD we have to per- 
form the microscopic large-n limit as mentioned 
above. It reduces to evaluating the rescaled (or 
unfolded) correlators 
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with a similar rescaling of the masses. Here, we 
have switched from positive to real Dirac opera- 
tor eigenvalues. It is easy to see that Eqs. (||) and 
([3]) result in the following expression for the mi- 
croscopic fc-point correlation function with (3Nf 
masses: 
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Here, pg' u '^ is the massless correlation func- 
tion which is entirely known Jl6[ in terms of de- 
terminants (/3 = 2) or quaternion determinants 
(/3 = 1,4). Since the universality of massless cor- 
relation functions has already been firmly estab- 
lished, the universality of the massive ones au- 
tomatically follows. Alternative representations 
of massive correlation functions were derived in 
Ref. [0 using Gaussian ensembles. There, the 
mass degeneracy for (3 = 4 is partially lifted to be 
two-fold. 

In the simplest situation of the spectral den- 
sity with (3 degenerate massive fermions, Eq. (0) 
reduces to 



P ( t' 0) (O 



P S (£,Wcc 



(8) 



The full density with (3 dynamical flavors is 
thus given by the quenched density p^'"'^ and 
the mass dependent correction term expressed 
through the connected part of the massless two- 
point correlation function pg"(£, ip)conn- 



The application of Eq. (||) to the symmetry 
class (3=1 was discussed in Rcf. and we will 
not consider it in what follows. However, it is 
instructive to consider the simplest example, the 
symmetry class (3 = 2. The connected part of the 
two-point correlation function is proportional to 
the square of the unitary kernel 

v ft \ - nJcx+i{r{)J a {i) fn s 

K cx{$,V) = ^772 T\ • ( 9 

2(r -?r) 

Combining Eqs. (|§|) and ([)[), we easily arrive at 
the result of Ref. jlO| for the microscopic density 
with two degenerate flavors, 



(10) 



Interestingly, the result of Ref. for non- 

degenerate masses p\ and pi can be put into the 
same form 



pf Va \i) 
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of the quenched density plus a mass-dependent 
correction. 

At (3 = 4, the microscopic density for four de- 
generate massive fermions is given by Eq. (^) 
with the massless microscopic density [flTj 
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and the connected part of the massless two-point 
correlation function || 

pf UA) (d,V)conn (12) 



where 



2? 



dtK 2v {t,2n) 



■2<l 



dtK 2v (2£,t). 



Performing the analytic continuation from J- 
Bessel to 7-Bessel functions completes the solu- 
tion of Eq. (|). 

In Fig. 1 the microscopic massive density 
pf vA) described by Eqs. (§), (0) and © is 
plotted for v = versus the lattice data of Ref. 
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Figure 1 . The microscopic density ps (£) plotted 
against lattice data for different values of \x. 



P3| with gauge group SU(2) in the fundamen- 
tal representation. Because of using staggered 
fcrmions symmetry class (3 = 4 applies. A reason- 
able agreement between our parameter-free theo- 
retical prediction and the lattice data is observed. 
The chiral condensate has been obtained from the 
Banks-Casher relation. A fit to the best value of 
E could improve the systematic shift for higher 
values of £ due to finite-size effect and statistics 
as is discussed in Ref. Jl3l . 
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